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1.  INTRODl TTION 


W hen  performing  computer  simukition>  of  continuous  and  discrete  dynamic  systems,  large 
amounts  of  output  data  are  often  gener.ited  and  presented  as  tables  or  plots  of  variables  as 
functions  of  time.  Since  much  of  the  program  writing  is  concerned  with  outputting  the  data, 
simulation  languages,  with  their  generalised  output  facilities,  olfer  many  advantages  to  the  pro- 
grammer. both  in  convenience  and  economies  of  elfort  lli>we\er.  an  ideal  form  of  output  is  not 
always  possible  because  of  limitations  in  the  variety  of  output  presentation. 

In  developing  a modilied  version  of  the  ‘block  iiriented'  simulation  language  CSMP-IO.' 
called  CSMP-IO(.\Rl.).-  one  of  the  major  changes  made  has  been  to  divide  the  language  into 
a modelling  program  named  HOMMP  (Block  Oriented  Mathematical  Modelling  Program) 
and  an  output  program  named  I R.ANS  ( I ranslation).  The  modelling  program  is  used  to 
perform  the  simulation  and  store  the  I'utput.  at  the  completion  of  each  time  step,  in  binary 
form  on  a specified  storage  device.  The  output  program  is  used  to  print  and  plot  the  character 
conversion  of  the  binary  tile  in  either  tabular  or  grapliic.il  form.  In  C'S.MF’-IO.  which  was 
developed  from  C‘SMP-9.-'  a signilicant  amount  id' core  storage  is  used  to  code  the  output  section 
of  the  program.  By  dividing  the  language  into  two  parts,  the  amount  of  core  storage  available 
for  each  process  is  significantly  increased.  This  allowed  more  facilities  to  be  included  in  both 
the  modelling  and  output  section  of  the  language  as  well  as  increasing  the  si/e  of  the  problem 
that  can  be  handled.  Because  the  'raw'  output  data  are  stored  in  a tile  in  binary  form  by  the 
modelling  program,  various  forms  of  output  of  the  same  dat.i  may  be  obtained  without  repeat- 
ing the  simulation.  A sav  ing  of  storage  is  made  in  the  binary  lile  by  storing  an  output  value  only 
if  it  changes  by  more  than  a prescribed  percentage  from  the  value  last  stored. 

In  Reference  2.  which  fully  describes  the  modelling  program,  only  a brief  account  of  running 
the  output  program  is  given;  a full  description  is  given  here.  Although  the  output  program  has 
been  written  primarily  to  be  used  in  conjunction  with  the  modelling  program,  the  output 
program  may  also  be  used  independently  as  a general  purpose  output  program,  provided 
data  are  supplied  in  the  appropriate  format.  ,A  with  BOM  Ml’,  the  program  TR  ANS  is 
designed  to  be  run  interactively  from  a leletvpe. 

Ihe  output  program  is  written  mainly  in  I ORfR.AN  l\  for  the  I’1)1’-I0  I ORTR.AN 
compiler  ‘I  40  \ ersion-27‘  running  under  the  operating  system  ‘TORO  fS'.  The  program  provides 
Teletype  and  line  printer  output  in  either  tabular  or  graphical  form  and  incremental  plotter 
output  in  the  form  of  'strip'  plots  (of  the  form  produced  by  a multi-channel  chart  recorder)  and 
'overlay'  plots  (a  single  graph  for  up  to  eight  output  variables). 

The  output  tile  obtained  by  running  BOMMP,  which  becomes  the  input  file  for  TR.ANS, 
includes  the  time  limits  and  interval,  a title,  and  labels  identil'y ing  all  the  variables.  Subject  to 
some  restrictions,  these  may  be  redefined.  I or  graphical  output,  convenient  scale  limits  are  cal- 
culated by  FRANS;  however,  these  limits  may  be  defined  by  the  user. 

The  basic  operations  of  the  output  program  necessary  to  produce  Teletype  output  are  given 
in  Section  2.  and  in  Section  ^ arc  described  the  operations  wliich  enable  the  user  to 
(u)  obtain  line  printer  and  incremental  plotter  outjnit. 

(/>)  redefine  the  time  parameters  and  labelling  information,  and 
(i  ) define  the  scale  limits  used  for  graphical  oulput 
With  the  aid  of  a complete  evample.  Section  4 shows  tiow  the  output  program  may  be  used 
without  running  the  modelling  program,  and  m Section  5 a description  is  given  of  the  structure 
of  the  computer  program  and  the  input  lile. 

2.  BASIC  OPERATIONS  OK  Ol  I PI  T PROCiRAM 


This  section  explains  how.  using  basic  operations,  the  output  program  may  be  run  to 
produce  Teletype  output.  Much  of  this  section  is  contained  in  Reference  2.  but  is  included  here 


for  completeness  of  presentation.  In  the  block  oriented  language  CSMP-I0(AR1-).  the  output 
variables  are  referred  to  as  blocks,  and  each  block  is  identified  by  a number.  The  file  containing 
the  output  values  of  these  blocks  is  therefore  referred  to  as  the  block  output  file  to  distinguish 
it  from  the  model  output  file,  which  is  used  by  the  modelling  program  to  store  the  various  state- 
ments specifying  the  current  model.-  The  output  (stored  in  the  block  output  file  named 
OUTPI  .DAT)  produced  by  the  non-linear  spring  problem  of  Reference  2 is  used  belowi  to 
illustrate  the  operations  of  TRANS.  Errors  and  their  diagnostic  messages  are  presented  in 
Appendi.x  A. 

2.1  Setting  up  the  Program 

For  the  PDP-IO  computer  the  relocatable  binary  file  TRANS. REL  is  first  loaded  and  a 
core  image  (TRANS. SAV)  is  saved.  E.xecution  of  the  program  is  then  initiated  by  the  user  and 
the  filename  (without  e.xtension)  of  the  block  output  tile  produced  by  running  BOMMP  is  typed 
by  the  user.  The  filename,  which  is  OUTPT  for  the  example  shown,  has  a maximum  of  five 
characters.  The  output  program  then  types  the  title,  date  and  time  of  creation  of  the  block 
output  file,  integration  interval,  run  CPU  time,  and  time  parameters,  followed  by  an  which 
signifies  the  command  mode  (see  below  ).  .^n  example  of  how  to  set  up  the  program  in  the  above 
manner  is  shown  below  :t 

Rl  TRANS 


IP  PILES  AMF.  OUTPT  lixtension  name ‘OAT‘ is  assumed 

SOS-I.ISEAR  SPRISO  HITII  DASH  POP 

IP  PILE  RECORPEI)  0\  :i-4PR-76  AP  lf>  .C’ 

l\TEO.\  IM  l.mHtE  DP  Rl  S ( PI  lIMP  1 SPC 

TIME  PRO  M OOODDE  Dl  PO  LODDDL  0/  l\  SIIP\  OI  I ihmhiL  do 


.Although  the  filename  extension  is  not  spesified  b\  the  user,  ihe  extension  name  ‘OAT’  is 
assumed.  If,  in  the  above  example,  a tile  named  Of  Iff  I>\1  is  not  fi'und  on  the  disk,  the 
message 


OI  TPT.DAT  S(yP  0\  DISK 
IP  IILES  IMP 

is  typed  by  the  program  so  that  a new  filename  may  be  specified. 

As  with  the  modelling  program,  the  output  program  is  controlled  by  commands,  which  are 
typed  by  the  user  whenever  an  is  typed  by  the  pri'gr.ini  ti  e whenever  control  is  at  the  com- 
mand mode).  The  commands  may  be  shortened  to  three  characters  .ind  are  listed  with  brief 
descriptions  in  Table  I 

By  using  the  appropriate  comni.inds,  one  or  more  p.iriicular  ty|X's  of  output  may  now  be 
specified.  To  process  the  input  tile,  producing  tabular  or  graphical  output,  the  ‘(iOl  ' command 
should  follow.  Control  is  then  returned  to  the  comiiiaiul  mode  Repetition  of  the  appropriate 
command  for  a particular  type  of  output  before  the  '(iOl  ' command,  deletes  the  effect  of  the 


t Teletype  messages  typed  by  the  user  are  shown  in  bold  upper  case  messages  type  by  the 
computer  program  are  shown  in  ifalic  uppercase  Comments  on  Teletype  messages  arc  shown 
in  upper  and  lower  case  print  alongstde  the  appropriate  message 


r 


TABI  F.  I 

Output  Prot>rani  Comniands 
All  commands  may  be  shortened  to  three  characters 


Comm.and 

Elfcct 

Section(s) 

for 

Reference 

FXIT 

Returns  control  to  the  monitor 

2.1 

goe: 

Processes  the  input  file  to  obtain  tabular  or  graphical  output 
for  the  commands  PRC Ol.UMN'.  PRPl.Or.  PLSTRIP'. 
and  PI  OVE  RLAY' 

2.1 

l.ABET. 

Eiiiables  modification  and  adililion  to  labelling  information 

3.2 

Pl.OVI  Rl  .S'! 

Specifies  'overlay'  plotting  on  the  incremental  plotter 

3.5 

PLSTRIP 

Specifies  'strip'  plotting  on  the  incremental  plotter 

3.5 

PRCOl.UMN 

Specifics  tabular  output  on  the  line  printer  or  releiype  in  the 
form  of  printed  columns 

2.2.  3.4 

PR PLOT 

Specifies  graphical  output  on  the  line  printer  or  Teletype 

2.2.  3.4 

RUN 

Fquivalent  to  the  commands  (in  order)  'PRC'Ol.UMN', 
'PRPLOT,  'PI  SI  RIP'.  and  'GOF';  if  any  form  of  output 
is  not  required,  a 'carriage-return'  is  typed  in  place  of  the 
block  numbers 

3.6 

.SC'ALT 

E nabics  tlie  user  to  specify  scale  limits  used  for  graphical 
output 

3.3 

time: 

Redefines  time  parameters 

3.1 

previous  use  ol’  that  command;  repetition  alter  the  ■CK)I.'  command,  results  in  additional 
output.  Control  is  returned  ti>  the  monitor  h>  the  Mv.XIT'  command. 

2.2  Teletype  Output 

The  commands  ‘PRCOl.UMN'  and  PRPLOT'  arc  used  to  obtain  tabular  and  graph/cal 
output  respectively  on  either  the  line  printer  or  lelelype.  On  the  graphical  output,  line  spaces 
correspond  to  equal  increments  in  time  (or  dependent  variable),  and  column  spaces  correspond 
to  increments  in  the  output  variables,  l or  each  of  these  comniands.  the  required  block  numbers, 
which  are  separated  by  any  non-numeric  character,  arc  first  listed.  They  may  be  typed  on  more 
than  one  line,  each  line  consisting  of  a maximum  of  72  characters,  and  the  list  is  terminated  by 
two  consecutive  ‘carriage-returns'.  If  all  the  block  numbers  arc  required,  the  user  may  type  A 
in  place  of  the  individual  block  numbers,  in  which  case  a single  ‘carriage-return’  terminates  the 
list.  When  specifying  blocks  by  their  individual  block  numbers,  a maximum  of  a hundred  bliKks 
may  be  processed  each  time  the  ‘GOT’  command  is  used;  however,  all  blocks  arc  processed  when 
the  user  types  A instead,  even  if  there  are  more  than  a hundred  blocks.  Following  the  Teletype 
message 

/.s  o r w n y ri-qrp  . 

the  user  types  \ if  the  output  is  to  be  printed  directly  on  the  Teletype,  or  any  other  character, 
including  a ‘carriage-return',  if  the  output  is  to  be  stored  on  a disk  file  so  that  the  resulting  file 
can  subsequently  be  printed  on  the  line  printer  (see  Section  .7.4).  This  latter  lorm  ol  output 
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contains  more  inl'ormalion  than  the  Teletype  output.  On  completing  the  Teletype  responses  to 
either  the  ‘PRCOl. LIVIN'  or  'PRPl.OT'  commands,  the  program  returns  control  to  the  command 
mode.  The  'GOf;'  command  is  then  used  to  obtain  the  Teletype  output.  It  should  be  noted  that 
both  tabular  and  graphical  output  cannot  be  obtained  on  the  Teletype  by  using  a single  ‘GOE’ 
command  following  the  'PRCOLUMN'  and  'PRPLOT'  commands. 

All  values  printed  on  both  forms  of  Teletype  output  are  shown  to  three  significant  figures. 
I or  the  Teletype  tabular  output,  up  to  six  variables,  in  addition  to  the  time,  are  printed  in 
columns  on  the  same  table.  If  more  than  six  are  to  be  printed,  six  at  a time  (or  the  remainder  if 
the  last  group  is  less  than  six)  arc  printed  in  separate  tables.  For  the  Teletype  graphical  output, 
up  to  live  variables,  represented  by  the  symbols  ' t '•■,  ‘ a and  ' are  displayed  on  the 

same  graph  with  a width  of  50  character  spaces  per  line  (i.e.  51  distinct  points).  Where  two  or  more 
symbols  coincide,  the  one  appearing  last  in  the  above  list  is  shown.  If  more  than  live  variables 
are  to  be  plotted,  live  at  a time  (or  the  remainder  if  the  last  group  is  less  than  live)  are  plotted  on 
separate  graphs.  Convenient  plotting  scale  limits  are  calculated  by  the  output  program  and  are 
shown  at  the  head  of  each  Teletype  plot  (sec  Section  for  specification  of  limits  by  the  user). 
The  output  is  printed  or  plotted  over  the  full  time  interval  for  which  data  arc  stored  in  the 
output  tile,  unless  the  'TIME'  command  is  used  to  restrict  the  time  interval  (see  Section  3.1). 

The  following  example  shows  how  both  tabular  and  graphical  Teletype  output  may  be 
obtained  for  the  non-linear  spring  problem: 

*PRC 

PRlXriM,  /.V  COLL  .UXS  : 

BLkS 

A l ist  terminated  by  a single  'carriage-return' 

IS  O P TO  TTY  REQRD  : \ 

*c;oE 

**  RL  .\\l\0  ** 

(See  Section  I of  Appendix  B for  Teletype  output) 

*PRP 

PRISTTR  PLOTS  : 

BLkS 

9.  4 I ist  terminated  by  two  consecutive  'carriage-returns' 

IS  O P TO  TTY  RTQRO  : V 

*c;oE 

•*  RL  SMSCi  •* 

( Sec  Section  2 of  Appendix  B for  Teletype  output) 

•EXI 


A .V/)  OP  LMX'l  nos 

CPL'  TIME  2. 27:  ELAPSE O TIME:  V.V.V2 

EXIT 


4 


3.  FURTHER  OPERATIONS  OF  OUTPl  T PROGRAM 


The  basic  operations  necessary  to  produce  Teletype  output  only  were  described  in  the 
previous  section.  This  section  describes  operations  enabling  the  user  to 
(a)  obtain  line  printer  and  incremental  plotter  output. 

(/))  redcline  the  time  parameters  and  labelling  information,  and 
(c)  define  the  scales  used  for  graphical  output. 

The  non-linear  spring  problem  output  is  again  used  as  an  example,  but  the  output  file  OUTPT.DAT 
is  obtained  by  running  BOMMP  with  an  output  interval  of  0.1  (previously  1.0);  this  allows 
smoother  curves  for  output  on  the  incremental  plotter.  Various  files  are  created  by  the  output 
program  to  store  labelling,  scaling,  or  output  iaformation.  The  labelling  and  scaling  information 
may  be  read  in  when  rc-runnmg  the  program,  and  the  output  information  enables  printed  and 
plotted  output  to  be  obtained  appropriately  on  the  line  printer  or  incremental  plotter.  The 
labelling  and  scaling  tiles  arc  named  TRANS. LBL  and  TR.ANS.SCL  respectively.  The  files 
containing  output  information  are  of  the  form  “filename". “extension",  where  “filename"  is 
the  filename  (without  extension)  of  the  block  output  file  produced  by  running  BOMMP  (i.e. 
OUTPT  for  the  present  example),  and  “extension"  is  set  by  TRANS  for  each  form  of  output. 
Table  2 provides  a summary  of  the  disk  files  created  by  TRANS  on  using  the  appropriate 
commands. 

TABl.E  2 

Disk  Files  Associated  with  Commands 
Files  containing  output  information  for  subsequent  printing  or  plotting 
have  the  filename  (without  extension)  of  the  block  output  tile  produced 
by  running  BOMMP  (i.e.  OUTPT  for  the  present  example) 


Command 

Name  of  File 

Logical 

Unit  Number 

1 ABFF 

TR  ANS  FBI. 

PLOVFRLAY.  PFSTRIP 

OL  rrn  pi.T 

9 

PRCOLUMN 

OUTPT  ( OF 

1 1 

PRPLOT 

OUTPT. PPL 

10 

SCALE 

TRANS. SCI. 

7 

3.1  Time  Parameters 

Preceding  the  output  values  stored  in  the  block  output  tile  produced  by  running  BOMMP 
are  the  lower  and  upper  time  values  and  the  output  time  interval  (sec  Section  5.1 ).  The  number  of 
time  values  for  which  printed  or  plotted  output  is  required  may  be  reduced  by  redefining  these 
parameters  using  the  ‘TIMF‘  command.  For  the  incremental  plotter  though,  the  ‘TIME'  com- 
mand has  no  elfcct  on  the  interval  used  for  plotting.  When  using  the  'TIME'  command,  times 
are  alwaxs  expressed  in  the  units  used  in  the  block  output  file.  For  the  present  example,  the 
U'wer  and  upper  time  values  of  0 and  20  respectively  are  unchanged,  and  the  output  interval  of 
0 I IS  altered  to  0.5  as  follows: 

*TIM 

r/t/A  l'  \R4\IS:  t.onUR.  IPPER.  ISTERPM.  (».  20.  0.5 
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.^.2  l.atK'lling 


Various  parameters  ami  labelling  information  are  stored  at  the  beginning  of  the  bIcKk  output 
file  produced  by  running  BOMMP  (see  Section  5.1).  On  reading  in  this  file,  TRANS  stores  the 
labelling  information  so  that  the  various  forms  of  line  printer  and  incremental  plotter  output 
produced  by  TRANS  may  be  .uitably  labelled.  The  labelling  information  includes  the  modelling 
program  title  of  one  line  of  sixty  characters  and  a list  of  block  labels  consisting  of  one  line  of  ten 
characters  for  each  block.  The  first  of  these  block  labels  is  the  independent  variable  label  (referred 
to  here  as  X label)  of ‘TIME’.  A block  number  of  —1  is  used  to  identify  this  particular  label. t 
Subsequent  block  labels  (referred  to  here  as  Y labels)  are  used  to  identify  the  block  variables 
(i.e.  dependent  variables).  If  no  label  is  specilied  when  running  BOMMP.  then  by  default, 
TRANS  provides  a label  of  the  form  ‘Bl.K  n’,  where  n is  the  appropriate  block  number.  For 
graphical  output,  blocks  are  identified  by  their  label,  rather  than  by  their  number  as  is  the  case 
with  Teletype  output.  Both  labels  and  block  numbers  are  shown  for  line  printer  tabular  output. 

The  ‘LABEL’  command  enables  the  title  and  labels  to  be  redefined  and  allows  the  speci- 
fication of  an  additional  line  for  each.  The  title  and  labels,  whether  modified  or  not,  are  stored  in 
a lile  named  TRANS. LBL  for  future  use.  On  using  the  ‘LABEL’  command,  the  labelling  infor- 
mation required  for  a particular  graphical  or  tabular  output  may  then  be  read  in  from  this  file. 
Labels  may  also  be  defined  for  blocks  whose  output  values  are  not  read  in  from  the  block  output 
file;  this  enables  the  user  to  set  up  a general  labelling  lile  for  use  with  a number  of  different  block 
output  files  w hich  each  have  a number  of  variables  in  common. 

Each  time  the  ‘LABEL'  command  is  used,  the  output  program  processes  three  separate 
stages.  In  order  of  execution,  these  are  termed  the  (u)  ‘source'.  (A)  'listing',  and  (<•)  'miHUftcaiion' 
stages.  The  start  of  each  of  these  stages  corresponds  to  the  following  Teletype  messages  typed 
by  TRANS 

(a)  IS  LABELUSG  TO  BE  READ  FROM  DSK  : 

(b)  IS  TTY  LISTISG  REQRD  : 

(c)  ARE  MODIFICATIOSS  REQRD  : 

The  user  answers  each  message  appropriately  by  typing  Y for  ‘yes'  or  any  other  character  for  ‘no’. 

For  the  ‘source’  stage,  if  the  user  responds  by  typing  the  labelling  information  is  read  in 
from  a file  named  TRANS. LBL;  otherwise,  the  labelling  information  already  stored  by  TRANS 
is  used. 

If  'listing'  is  required,  the  program  first  types  the  title  followed  by  the  message 
BLKS  FOR  LISTIS'G 

The  user  then  specifies  the  appropriate  block  numbers,  either  individually  on  a single  line 
consisting  of  a maximum  of  72  characters  or  by  typing  A if  all  the  block  numbers  arc  required. 
In  each  case,  a single  ‘carriage-return’  terminates  the  list.  If  individual  block  numbers  are  specified, 
following  the  listing  of  the  corresponding  labels,  the  above  Teletype  message  is  repeated  so  that 
further  block  numbers  may  be  specified.  This  priKess  is  repeated  until  the  user  types  a ’carriage- 
return’  only  in  place  of  a list  of  block  numbers.  If  all  the  block  numK*rs  arc  specified  (i.c.  by- 
typing  A),  the  listing  section  is  automatically  terminated  following  the  complete  listing  of  the 
labels. 

If ‘modifications’  are  required,  then  for  each  line  of  labelling  information  to  be  modified, 
the  user  types  the  block  number  followed  by  the  line  number  and  labelling  information  text.  For 
the  labels,  the  block  number  is  set  equal  to  the  number  used  to  identify  the  corresponding  output 
value,  but  for  the  title,  the  block  number  is  set  equal  to  zero.  I or  both  the  title  and  labels,  the 
line  number  is  set  equal  to  one  if  the  labelling  information  of  one  line  read  in  from  the  block 


■f  In  the  modelling  program,  block  number  I is  reserved  for  the  time  or  independent  variable. 
Because  the  output  program  may  be  used  independently  of  the  modelling  program  (see 
Section  4),  TRANS  identifies  block  number  I as  the  time  or  X value. 
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output  flic  is  to  be  modified,  and  is  set  equal  to  two  if  the  additional  line  is  to  be  specified  ot 
modified.  When  the  modifications  are  completed,  control  is  returned  to  the  command  mode  by 
typing  an  additional  ‘carriage-return'. 

The  two  examples  given  below  show  how  (I)  the  title  and  labels  stored  in  the  block  output 
file  OUTPT.DAT  arc  first  listed  on  the  Teletype  and  then  modified,  and  (2)  by  repeating  the 
'LABliL'  command  during  cither  the  same  run  or  a subsequent  one,  the  modified  title  and  labels 
arc  read  in  from  the  file  TRANS. LBL  and  listed  on  the  Teletype.  In  Example  2,  the  same  listing 
would  be  obtained  without  reading  the  labelling  information  from  TRANS. LBL  (i.e.  the  user 
types  N to  the  above  message  (a) ) if  the  ‘LABEL’  command  was  repeated  during  the  same  run. 


Example  / 


M.AB 

BLK  SO.  -I  DESOTES  ISDEP  V.ARl.ABLE 
IS  LA  BELL!  SO  TO  BE  READ  EROM  DSK  : N 
IS  TTY  I.ISTI.SG  REQRD  V 
TITLE 

NON-LINEAR  SPRING  WITH  D.ASHPOT 

Additional  line  is  blank 


By  typing  A,  ‘listing'  stage  automatically 
terminated  following  complete  listing 

Additional  line  for  each  label  is  blank 

4 ) DBLE  DOT 
4S  y DOT 


BLKS  FOR  LISTISC 
A 

BLK  LABEL 
-/  TIME 

9 Y 


10  BLK  a 10 


ARE  MODIEICATIOSS  REQRD  : \ 

BLK  (0  FOR  TITLE).  LISE.  TEXT 

0,  1,  NON-LINEAR  SPRING  WITH  HASHPOT 

0.  2.  DATA  (STORED  IN  OL  TPT.DAT)  PRODUCED  BY  RUNNING  BOMMP 

-1.  2,  SEC 

4.  2,  FT/SEC/SEC 

9.  2,  FT 

4«.  2.  FT /SEC 

Additional  'carriage-return'  terminates  the  list 
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F.xamplf  - 


*I.AB 

lil.K  \().  I PI-SOrFS  IS'DFP  VARIABLE 
IS  I.ABEI.LISli  TO  BE  REAP  FROM  PSK  \ 
IS  IT)  I.ISTISG  REQRP  : \ 


TIFLE 


SOS  I.ISEAR  SPRIXO  H U II  PASIIPOT 
DATA  (SrOREP  IS  OVTPT.PAT)  PROPLCEP  BY  RL  SMSG  BOMMP 


BITS  FOR  I ISTISG 

-I,  4.  •).  48  Single  'earriage-rcturn'  lerminates  the  list 

BIT  LABEL 

-I  TIME 

SEC 

4 ) PBLE  POT 
FT  SEC  SEC 
V ) 

FT 

4S  ) POT 

FT  SEC 

BITS  FOR  I.ISTISG 

‘Carriage-return"  at  beginning  of  line  lermin- 
ARE  .MOPIFICATIOSS  REQRP  : N ates ‘listing"  stage 


The  hie  TRANS. I BL  ma>  be  created  or  niodilied  without  running  TRANS.  A listing  of 
TR.ANS.l  BL  created  in  Txaniple  2 above  is  given  below  to  illustrate  the  exact  form  of  the  hie; 
the  format  statement  is  K)RM,A1(I4,  IX.  12A3)  for  each  line  of  the  title,  and  FORMAT(I4, 
IX.  4A5)  for  each  label: 

0 NON-1  INE i.AR  SPRING  WITH  DASHPOT 

0 DATA  (STOREiD  IN  OUTI’TDAT)  PRODUCE  Ei  ED  RUNNING  BOMMP 

1 teme:  set 

4 \ ET 

4 V DBI.E;  dote  T SFC  SE  C" 

48  DOE  FT  SFC" 

10  BEK//  10 

.V.4  Scaling  for  Cirapliical  Output 

On  reading  in  the  block  output  hie  produced  by  running  BOMMP.  the  output  program 
determines  lower  and  upper  data  limits  for  the  independent  variable  (referred  to  as  X)  and  for 
each  of  the  block  variables.  Whenever  the  "TIME:"  command  is  used,  the  data  limits  arc  re- 
determined for  the  reduced  number  of  data  values.  When  graphical  output  is  required,  these 
limits  are  then  used  to  calculate  convenient  plotting  scale  limits.  The  ‘SCALE"  command  enables 
this  process  to  be  bv passed  and  the  scale  limits  to  be  provided  by  the  user  and  stored  in  a hie 
named  TR.ANS.SCE.  for  future  use.  Eiata  limits  determined  by  TRANS  on  reading  in  the  block 
output  tile  arc  not  stored  in  TRANS. SC"1..  On  subsequently  running  TRANS,  the  limits  required 


fora  parlkular  graphical  output  may  be  read  in  from  the  lile  TRANS. SCL  on  using  the  ‘SCALE’ 
command.  As  with  labelling  information,  lower  and  upper  limits  may  be  defined  for  blocks 
whose  output  values  are  not  always  read  in.  When  the  user  defines  the  scale  limits,  they  are  used 
for  graphical  output  w ithout  any  modification.  Scale  limits  defined  for  the  X scale  are  used  only 
for  incremental  plots,  whereas  graphical  output  on  the  Teletype  or  line  printer  uses  the  time 
limits  specified  in  the  input  file,  or  the  modified  limits  if  the  ‘TIME’  command  is  used. 

Use  of  the  ‘SC.ALE’  command  is  very  similar  to  that  of  the  ‘I.ABEL'  command,  and  the 
two  examples  below  correspond  very  closely  to  those  used  to  illustrate  the  ‘LABEL’  command. 
Ditferences  may  be  seen  by  comparing  the  examples.  F’or  the  ‘listing’  stage,  the  limits  are 
followed  by  ‘ -v  if  they  are  data  limits  determined  by  TR.ANS  on  reading  through  the  block 
output  file.  F or  modifications,  on  each  line  the  user  types  the  block  number  followed  by  the 
lower  and  upper  scale  limits. 


Example  / 


*SC  A 

BLK  SO.  -I  DESOTES  IS  DEE  \ .A  RI.-\BLE 


.ARE  PLOT  SCALE  LIMITS  TO  BE  RE.AD  EROM  DSK  : N | 

I 

I 

IS  TTY  I.ISTISC,  Ol  LIMITS  REQRO  : \ 


^ DESOTES  DATA  LIMITS  OSLY  {LSED  TO  DETERM  IS  E CLOT  SC. ALES) 


BLKS  TOR  LIST  ISO 
A 


BLK  LOWER 
— I O.dOOOE  01 
9 -I.OOOOE  01 
4 S.75I4E-00 

4S  -4.7SSSE  ■ 00 
10  I.OOOOE  02 


L PEER 

2. 0000 E 01  ■■■  ■■■■  ■■■' 
6.46S<7E  00  .v 
2. 0000 E 01 

>4.4473 E ■ 00  // 

4.2564E  01  -v 


ARE  MODITIC.ATIOSS  REQRD  : \ 

BLK.  LOW  ER.  L EPER 

10,  -60.  40 
4.  -10.  14 


Example  2 

*SC\ 

ELK  SO.  I DESOTES  IS  DTE  T.ARLABI  E 

ARE  PLOT  SCALE  LIMITS  TO  BE  READ  PROM  DSK  . > 

IS  TTY  I.ISTISC,  Ol  LIMITS  REQRD  N 

-•/  DESOTES  DATA  LIMITS  OSLY  {USED  TO  DETERMISE  PLOT  SCALES) 
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r 


BLKS  FOR  I.ISTIS'C 

4,  10,  48 


HI  K 

LOWER 

1 PPER 

4 

1.0000  E 01 

I.4000E  ■ 01 

10 

6.0000E  01 

4.mtOE  ■ 01 

4S 

4.7.5S,V/;'  00 

0.4473E  ■ 00  If  4 ^ 

RLFS  FOR  LISTIMi 

ARE  \fO/)lFlCAFIO\S  REQRD  . N 

♦ 


As  with  TRANS, LBL,  the  tile  TRANS. SCL  may  be  created  or  modilied  without  running 
TR  ANS.  A listing  of  TR  ANS. SCL  created  in  Example  2 above  is  given  below  to  illustrate  the 
exact  form  of  the  file;  the  format  statement  for  each  line  is  FORVIAT(14,  2(1  PE  1 2.4)  ).  and 
because  the  limits  for  only  two  blocks  were  defined,  the  file  contains  only  two  lines: 

4 l.0(X)0E  01  1,4000E  01 
10  - 6.0000E  ■ 01  4.0000E  01 

3.4  Line  Printer  Output 

The  commands  'PRCOLUMN'  and  'PRPLOT'  used  to  obtain  tabular  and  graphical  output 
respectively  have  been  described  in  Section  2.2.  However,  only  the  Teletype  form  of  output  was 
described  in  detail.  The  line  printer  form  of  output  obtained  from  the  disk  file  is  fully  described 
in  this  section.  The  tabular  and  graphical  outputs  for  the  line  printer  are  stored  by  the  output 
program  in  the  files  "filename". COL  and  "filename". PPL  respectively  ; both  forms  of  output 
may  therefore  be  obtained  on  using  a single  ‘GOE’  command  following  the  ’PRCOLUMN'  and 
■PRPLOT  commands.  On  the  line  printer,  all  values  printed  are  shown  to  five  significant  figures 
and  labels  are  used  to  identify  the  blocks  instead  of  the  block  numbers  for  graphical  output,  and 
both  labels  and  block  numbers  are  show  n for  tabular  output.  The  title,  date  and  time  of  creation 
of  the  block  output  lile.  integration  interval,  and  run  CPU  time  arc  also  printed.  For  the  tabular 
output,  up  to  ten  variables,  in  addition  to  the  time,  are  printed  in  columns  on  the  same  table, 
and  for  the  graphical  output,  the  graph  has  a width  of  100  character  spaces  per  line  (i.e.  101 
distinct  points)  and  symbols  are  overprinted  by  other  symbols  for  points  that  coincide. 

The  following  example  shows  how  the  tiles  OUTPT.COL  and  OUTPT.PPL  may  be  generated 
for  the  present  example;  the  redefined  time  parameters,  labelling  information,  and  scales  specified 
in  Sections  3.1  to  3.3  are  used; 


•PRC 

FRISTISO  /A  COLL  MSS  : 

m.KS 

A 

IS  OP  ro  TTY  REQRD  : N 

•PRP 

PRIS TER  PLOTS  : 
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BLKS 

A 

/S  O P ro  TTY  RPQRD  : ^ 
*GOK 

**  RVSMSC,  ** 


On  .subsequently  using  the  ‘EXIT'  coinmaml  in  the  above  example,  the  tiles  OUTPT.COI. 
and  OUTPT.PI’E  may  be  printed  b\  running  the  syst'-m  program  'PRINT',  and  are  shown  in 
Appendix  C.  Because  carriage  control  characters  are  Used  extensively  to  produce  the  graphical 
output,  the  file  OUTP1.PPL  should  be  printed  using  the  'P'  switch  (on  the  PDP-IO  computer) 
as  follows: 

.R  PRIM 


/7.-/(t*()lTPT.PPI./PS 

rnj:s  PRisTiin  req.  pages 

OLTPT  PPl.  ROI  2 

EMI 


C arriage  control  characters  are  used  minimally  to  produce  the  tabular  output  so  that  use  of 
the  'P'  switch  is  optional  when  printing  the  file  OUTPT.COI.. 

3.5  Incremental  Plotter  Output 

The  commands  ‘PESTRIP'  and  ‘PLOVT  RI.A^  " are  used  to  obtain  'strip'  plots  (of  the  form 
produced  by  a multi-channel  chairt  recorder)  and  'overlay'  plots  (;i  single  graph  lor  up  to  eight 
output  variables),  liach  form  of  output  Is  stored  in  the  tile  “filename  Pl  f.  Time,  or  the  inde- 
pendent variable,  is  shown  on  the  aibscissa  (.\  aixis).  and  the  output  variables  are  shown  on  the 
ordinate  (V  axis).  .As  with  line  printer  output,  the  title,  date  and  time  of  creation  of  the  block 
output  tile,  integration  interval,  and  run  C PU  time  are  shown  on  each  plot.  Labels  arc  used  to 
identify  the  blocks  instead  of  the  block  numbers. 

Tor  both  commands,  the  block  numbers  required  are  first  listed  in  the  same  way  as  for  the 
'PRC  C)l. LAIN"  and  'PRPI  Of  commands  (i  e.  the  user  types  A or  specifies  individual  Hock 
numbers).  Tor  the  'IT  ()\T  R1  AN  ' command  only,  the  block  numbers,  whose  lower  and  uppt. 
variable  limits  are  used  to  determine  the  'i'  phitting  scale  arc  then  listed  m the  same  way.  Finally, 
for  both  commands,  the  lengths  I'f  the  \ and  N axes  expressed  as  an  integer  number  of  inches 
are  then  spccilied;  the  length  v'f  the  N axis  must  not  exceed  8 in. 

I he  'strip'  plots  are  .luti'ni.itically  arranged  in  adiaccnt  strips  such  that  the  total  span  of  the 
Y axes  does  not  cxcceil  8 in.  further  sets  of  plots  are  createii  until  all  the  variables  arc  plotted. 
As  with  graphical  I'utput  on  the  leletyix’  or  line  printer,  cv’nvcnient  plotting  scale  limits  are 
calculated  bv  the  progr.im  for  each  output  variable,  unless  the  sc.ile  limits  are  specified  by  the 
user.  Tor  the  'overlay'  plots,  up  to  eight  v.iri.ibles  may  be  plotted  on  the  same  graph  using  a 
c vmmon  scale  If  there  are  more  than  eight,  only  the  tirsi  eight  are  plotted;  the  remainder  are 
ignored  anil  an  appropri.ite  error  message  is  typed  (see  .Appendix  .A).  LIsing  only  the  block 
numbers  specitied  for  determining  the  N plotting  scale,  the  smallest  lower  limit  and  largest  upper 
limit  are  determineil  It  both  limits  .ire  scale  limits  specitied  by  the  user,  they  are  used  without 
modilicalion  lf.it  lea-t  one  limit  is  .i  data  limit,  a convenient  plotting  scale  is  calculated  by  the 
program.  Including  the  title  .mil  l.ibels.  a plot  will  lit  within  the  lntcrnati<>nal  standard  .A4  paper 
si/c  for  an  .\  axis  of  length  5 in. 
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The  two  examples  given  below  show  how  for  the  present  example,  using  the  modifications  of 
Sections  3.1  to  3.3.  the  tile  OUTPT.PLT  may  be  obtained  for  (I)  'strip'  plots  and  (2)  'overlay' 
plots.  Two  separate  runs  of  TRANS  were  made;  however,  both  plots  may  be  obtained  on  the 
one  file  by  using  a 'COE'  command  after  each  plotting  command. 

Example  I 

*PLS 

STRIP  PLOTS  : 

BLES 

4,  48.  9,  10  Full  list  of  block,  numbers  typed  individually  in  place  of  A 

to  obtain  desired  order  of  plots 
LESGTH  OE  AXES  /.V  IXCHES:  X,  >'  5,  2 

*GOE 

♦♦  RU.WIXG** 

* 


Example  2 
*P1.0 

OTERI.A  ) PLOTS  : 

BLES 

9,  48,  4 

BLES  TO  OETERMIXE  > SCALE 

4 List  terminated  by  two  consecutive  'carriage-returns',  as 

in  previous  list 

LESGTH  OE  A XES  L\  IXCHES;  .\ , >'  5,  8 

*GOE 

**  RUXXIXG  *• 


* 


On  subsequently  using  the  'EXIT'  command  in  each  of  the  above  examples,  the  file 
OUTPT.Pl.T  may  be  plotted  (see  Figures  la  and  lb)  by  running  the  system  program  'PLOTQ' 
as  follows: 


.R  PI.OTQ 


/().37*Ol  TPT.Pl.TS 
TILES  PLOTTED  REQ. 
OLTPT  PIT  Rdl 
EXIT 


EEET  LIMITS 

2 XtS.Atr.  T 1(1X0').  ) ( ft  ) 


L 


w 


3.6  The  ‘Rl!N’  Command 

The  ‘RUN’  command  is  equivalent  to  the  individual  commands  'PRCOLUMN',  ‘PRPLOT’, 
■PLSTRIP'.  and  ‘GOP’  in  the  order  shown.  If  any  form  of  output  is  not  required,  a 'carriage- 
return'  is  t>ped  in  place  of  the  block  numbers.  The  use  of  the  'RUN'  command  is  demonstrated 
by  the  following  esample,  which  repeats  the  examples  of  Sections  .^.4  and  3.5  (except  for 
'overlay'  plotting): 

♦RUN 

PRI.STI.SG  l.\  COLL  . MSS  : 

Hl.KS 

A 

/.S'  O r TO  TTY  RKQRD  : N 

PRISTLR  PLOTS  : 

BLKS 

A 

IS  O P TO  TTY  REQRD  : N 
STRIP  PLOTS  : 

BLKS 

4,  4«.  9,  10 

LESGTH  OF  .4\ES  IS  ISCHES:  X.  X = 5.  2 
**  RUSSI.SG  ** 


* 

4.  I SK  OF  on  PI  T PROGRAM  WITHOUT  Rl  NMNG  THK 
MODELLING  PROGRAM 

A subprogram  package  consisting  of  the  FORTRAN  file  PARAGE. F4  and  MACRO-IO 
tile  CPU. MAC  is  available  to  enable  the  user  to  obtain  an  input  tile  for  TRANS  without  running 
BOMMP.  An  example  of  a main  program  that  generates  certain  functions  and  stores  the  results 
in  the  file  FUNCT.DAT  in  binary  form  for  input  to  TRANS  is  given  in  Appendix  D.  Comments 
in  the  listing  of  the  example  main  program  provide  a complete  description  on  how  to  use  this 
facility.  The  following  is  an  example  run  of  TRANS  that  uses  the  input  file  FUNCT.OAT 
created  on  running  the  above  program: 


.RU  TRANS 

/ P I ILFS. 4 ME  Fl'NC  T 

\ .4RIOL  S EL  SCTIOSS  OF  X 
IP  FILE  RECORDET)  OS  27-JIL-76  .47  N.o: 

I.XTEGS  1ST  D.OOOOE  01;  RL'S  CPC  TIME  1.15  SEC. 


TIME  FROM  O.OOOOE  01  TO  T IMHIOE  00  IS  STEPS  OF  4.(HmE  O.i 


LAB 


BLK  SO.  -I  Pf:.\OTt:S  l.\DEP  VARIABLE 
IS  LABELLISG  TO  BE  READ  FROM  DSK  : N 
IS  TTY  LISri.\G  REQRD  : Y 
TITLE 

V A RIO  VS  FUNCTIONS  OF  X 

BLKS  FOR  LISTING 
A 

BLK  LABEL 
-I  A 

1 EXP{2X) 

2 SQUARE 

3 SIN(S\) 

4 X*SL\{,'i  \ ) 

5 SIS3SX**2) 

6 SAM 

7 NORMAL 

ARE  MODIFICATIONS  REQRD  : V 

BLK  {()  FOR  TFII.E).  LINE.  TEXT 

2.  2,  WAVK 
6,  2,  TOOTH 
7,2.  riRVK. 

•SC  A 

BLK  NO.  I DE.SOTES  INDEP  VARIABLE 

.\RE  PLOT  SCALE  LIMITS  TO  BE  READ  FROM  DSK  : N 

IS  TTY  LISTING  OF  LIMITS  REQRD  : Y 

■V  M 24  DENOTES  D.4TA  LIMITS  ONLY  (USED  TO  DETERMINE  PLOT  SCALES) 
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BLKS  FOR  IJSTIS'G 
A 

BLK  LO  IFF  R UPPER 

- 1 O.mWE- 01  2mH)E  ■ IH)  # # # 

/ 0.(MHM)E-0l  5.459HE  ■ 01  ### 

2 _ I.0(HH)E  ■ 00  I.mME  ■.  (X)  # # # 

J -9.97JRE  01  <4.V777E~OI 

4 ~IJ75'7E  ■ (X)  1.757 IE  • (X)  # # 

5 ~9.WWE  01  9.9m)E~0l  77  77-77 

6 0.(XXX)E  0!  4.mX)E-0I  f7. 77  77 

7 0.(XXX)E-OI  3.<7m)E-0l  77  77  77 

.4 RE  MOniEIC.4TIO.SS  REQRl)  Y 

BLK.  LOWER.  UPPER 

1.  0.  20 

*PLS 

STRIP  PLOTS  : 

BLKS 

A 

LESGTH  OF  A \ES  IS  f.SCHES:  -V,  > 5,  I 

•GOE 

♦*  RU.S.SLSG  ** 

♦EXI 


E.SD  OF  E.SECUTIOS 

CPU  TIME:  X.60  ELAPSED  TIME:  4:49.52 

EXIT 

.R  PI.OTQ 


/4  ’-/’El  NC  T.PETS 

FILES  PLOTTED  REQ.  FEET  LIMITS 

l U.SCT  PLT  ROI  2 X{S.50’.  F 10.50').  Y(-5.32',  5") 

EXIT 


The  output  file  FUNCT.PLT  obtained  above  is  shown  in  Figure  2. 

In  the  above  example,  the  labelling  information  and  plotting  scales  were  modified  inter- 
actively while  running  TRANS.  This  may  be  inconvenient,  especially  if  frequent  use  of  the 
additKvnal  line  fi>r  the  title  or  labels  is  required.  To  avoid  this  inconvenience,  the  modifications 
may  be  elVccted  directly  in  the  main  program  by  opening  the  files  TRANS. FBL  and  TRANS. SCL, 
and  storing  the  modified  information  in  these  files  in  the  format  specified  in  Sections  3.2  and  3.3. 
On  running  TRANS,  the  labelling  information  and  plotting  scales  would  then  be  read  in  using 
the  T ABF  I,'  and  ‘SCAl.F'  commands  respectively. 
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5.  DFSrRlPTION  OF  Ol  TPl  I PROGRAM 

1 he  input  tile  ^truclllrc  and  computer  program  structure  are  hrlellv  described  below,  l urther 
information  may  be  obtained  by  referring  to  the  FORI  R.AN  listing. 

5.1  Input  File  Structure 

The  input  file  for  TR.ANS.  which  is  generally  the  block  output  (ilc  obtained  by  running 
BOMMI*.  is  a binary  tile  consisting  of  full  disk  blocks  of  I2S  words.  The  modelling  program 
arranges  the  data  into  records  of  126  data  words  with  two  additional  words  being  used  by  the 
operating  system  'FOROTS'  to  signify  the  beginning  and  end  of  each  record 

To  save  storage,  when  running  BO.M.MPan  output  \alue  is  only  stored  if  it  changes  by  more 
than  a prescribed  percentage  from  the  value  last  stored  at  .i  previous  time  value.  Ihe  data  are 
therefore  stored  in  the  form  of  pairs  consisting  of  a block  number  (stored  as  a real  variable 
rather  than  as  an  integer)  followed  by  the  corresponding  output  value  A block  nuntber  equal  to 
— I is  used  to  indicate  that  the  next  word  is  a new  time  or  independent  variable  value.  For  N 
blocks,  the  first  24  .^N  words  of  the  input  tile  are  used  to  stiirc  v.irious  parameters  and 

labelling  information  (sec  Appendix  Ei).  The  next  2 2N  words  consist  of  a block  number  of  — 1. 
the  lower  time  value,  followed  by  N block  number  and  block  value  p.iirs.  Subsequently,  a block 
number  of  I.  a time  value,  followed  by  a maximum  of  N block  number  and  value  pairs  are 
stored,  the  actual  number  depending  on  whether  the  block  values  have  changed  sufficiently  from 
their  value  last  stored  at  a previous  time  step.  Finally,  the  run  C PL'  time  is  stored. 

5.2  F’rogram  Structure  and  Storage  Requirement 

The  computer  program  consists  of  a small  main  program  (MAIN  ) used  to  open  the  arrays 
that  are  later  expanded,  a number  of  major  subroutines  and  various  service  subroutines  and 
functions  (see  Appendix  F).  When  loaded  on  the  PDP-10  computer  the  program  requires  I.E  K 
words  of  core  store  ( I K 1024).  On  running  the  program,  a number  of  arrays  arc  expanded  so 
that  an  additional  8N  words  are  required,  where  N is  the  number  of  blocks. 

6.  C ONCT.l  DING  RKMARKS 

Although  designed  to  be  used  as  the  output  section  of  the  'block  oriented'  simulation 
language  C.SMP-I()(.ARI.).  the  output  program  that  has  been  described  may  also  be  used  as  a 
general  purpose  output  program  by  using  a supplied  program  package.  .As  with  the  modelling 
program,  which  is  used  to  perform  the  model  simulation  in  C'SMP-10(.ARI  ) and  store  the  output 
in  a binary  tile,  the  output  program  is  written  mainly  in  FOR  FR  AN  l\'  for  a PDP-IO  computer 
and  is  controlled  interactively  from  a Feletype.  Provision  h.is  been  made  for  Teletype  and  line 
printer  output  in  either  tabular  or  graphical  fi>rm  and  incremental  plotter  output  in  the  form  of 
'strip'  or  'overlay'  plots.  Facilities  have  been  introduced  that  provide  a large  degree  of  flexibility 
in  allowing  the  user  to  redefine  the  time  parameters,  labelling  information,  and  plotting  scales. 
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APPKNDIX  A 

Output  Program  Diagnostic  Krror  Messages 

To  distinguish  the  following  error  messages  from  the  TOR  I RAN  compiler  diagnostics,  all 
the  messages  below  are  preceded  by  an  exclamation  mark  ( !).  An  integer  constant  is  represented 
by  n and  's.e.'  is  used  to  abbreviate  'self-explanatory'. 

CODE  PRIOR  TO  CALI.  TO  EXPAND  UNRECOGNIZED  obtain  a new  copy  of  the 
file  I R.ANS.S.AV;  if  error  persists,  keep  data  files  and  consult  originator 

COM M AS  D ERROR  invalid  program  command 

END  OE  IP  TILE  REACHED  BEEORE  UPPER  TIME  TALUE  s.e. 

EH.E  TRANS  I III.  {CONTAINING  LABELLING  INEOR.MATION)  NOT  ON  DSN — 
labelling  information  already  stored  by  TRANS  is  used 

TILE  TRANS. SCI.  (C'0.\r,-1/,\7.V(;  SCALE  LIMITS)  NOT  ON  DSE  data  limits  deter- 
mined by  TR.ANS  on  reading  in  the  block  output  file  are  used 

EIRST  too  BLOCKS  ONLY  HILL  BE  PROCESSED  when  specifying  blocks  by  their 
indisidual  block  numbers,  a maximum  of  a hundred  blocks  may  be  processed  each  time  the 
'GOIf  command  is  used  following  one  or  more  of  the  commands  'PRCOLUMN',  'PRPLOT' 
and  PLSTRIP' 

EIRST  (V  BI.KS  ONLY  HILL  BE  PROCESSED  a maximum  of  eight  blocks  may  be 
processed  for  'overlay'  plots  each  time  the  'GOF.'  command  is  used  following  the 
'PLOVERl.AV  ■ command 

INCOMPLETE  TILE  for  incorrect  or  incomplete  format  for  parameters  and  labelling 
information  at  beginning  of  input  tile  (see  Appendix  K);  may  occur  when  creating  the  input 
file  without  running  BOMMP  (see  Section  4) 

LENGTH  OE  .\  AMS  MUST  BE  A I E INTEGER  s.e. 

LENGTH  OE  I AMS  MUST  BE  A YE  INTEGER  LESS  THAN  9 —s.e. 

.'  I.OH  ER  LIMIT  MUST  BE  LESS  THAN  UPPER  LIMIT  when  specifying  scale  limits 

.'  NOT  ENOUGH  CORE  TO  SET  UP  ARRA  YS  when  expanding  arrays  in  'EXPAND' 

.'  ONE  )■  / EORM  OE  O P TO  TTY  PER  RUN  both  tabular  and  graphical  output  cannot 
be  obtained  on  the  Telety  pe  by  using  a single  'GOE'  command 

! O P EOR  BI.K  n NOT  STORED  IN  / P ITLE  when  specifying  blocks  to  be  output 

.'  SYNTA.y  ERROR  typing  error  or  inappropriate  response  to  Teletype  message 


IS 


AIM’HNDIX  B 
Teletype  Output 

( 1 ) Tabular  Output 


TIME 

e.eec-ei 

1. aeE»M 

2. eec*ee 

3. ME»0e 

4. eeE*ee 
s.eeE«ee 
6.eeE*ee 
7.00E*ee 
8.e0E«ee 

9.00E*00 
I .00C*0t 
I . I0E*0I 
! .20E*0I 
I .3eE*0l 
I .40E*0I 
I .S0E«0l 
I . 60E*0I 
1 .70E*01 
I .80E*01 
I .90E«0I 
2.00E*0I 


BLK*  9 

-I .eeE»0i 
•3.27E«0e 
4.5eE*0e 
6.06E«00 
I .86E*ee 
-2.02E«0« 
-3.68E*0e 
-2.69E*0e 
-7.47E-0I 
8.00E-0I 
I . S8E*0e 
t . 60F*0e 
I . i0E«0e 
4. I 7E-0I 
-I.75F-0I 
-5. 1 6E-01 
-5.80E-0I 
-4.37E-0I 
-2.02E-01 
I .96E-02 
I .62E-0I 


BLK*  4 
2.00E«0I 
-l.42E«00 
•6. 41 E«00 
-6.39E*0B 
I . I0E*00 
2.0IE*00 
2.90E* 00 
8.9SE-0I 
-4.38E-0I 
-7.97E-B! 
-7.82E-01 
-5.23C-0I 

- I .83E-0I 
I .04E-0I 
2. 63F-0I 
2.85E-0t 
2.09E-0I 
9. 1 6E-02 

- 1 . 64E-02 
-8.35E-02 

- I .02E-0I 


BLK#  4B 
0.00E-0I 
9. 36E*00 
S. 52E*00 
•2. 44E*00 
•4.62E*00 
-2.98E*00 
-2.03E-et 
I .84E*00 
I .84E*00 
I . I9E«00 
3. 74E-0I 
•2.92E-0I 
-6.46E-ei 
-6.76E-0I 
-4. 8 IE- 01 
- 1 .97E-01 
S. 66E-02 
2.08E-01 
2. 43E-0I 
I .89E-01 
9.25E-02 


BLK#  IB 
-l.0eE«B2 

-I.I6E*0I 
2. I0E«0I 
3.  68E«BI 
3.7IE*0e 
-4. I0E«00 
-I .4IE*0I 
-8« ISE*00 
- I.49E*00 
l.60E*00 
3« 1 cE^BB 
3.20E*00 
2.20E*00 
8.33E-0I 
-3. 51 E-01 
- t .03E*00 
- I . I 6E*00 
-8.74E-0I 
-4.05E-01 
3.92E-02 
3.24E-0I 


(2)  (iraphical  Output 


BLK#  9 I -I.00E+0I  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  l.0eE*0t 

BLK#  4 I -2.00E*0I  •••••*•••••••*••••**  2.00E»0I 


TIME  I 


0.00E-01  * 
l.00E«00  I 
2.00E«00  I 
3.00E«00  I 
4.00E«00  I 
S.00E*00  I 
6.00E*00  I 
7.00E«00  I 
e.00E»00  1 
9.0eE*00  1 
l.00E*0l  I 
I.I0E»0I  I 
l.20E*0l  I 
l.30E«0l  1 
l.40E«0l  1 
I.S0E*0l  I 
l.6ec«0l  I 
l.70E*0l  I 
l.80E*0l  I 
l.90E*0l  I 
2.0eE»0l  I 


APPKNDIX  C 
l.inc  Printer  Output 


(1)  Tabular  Output  {File  OLTPT.COL  printed  u.sinf^  'P'  switeh) 


NON-LI’^EaH  sowing  with  DAShhdT 
data  (STORES  IN  OJTPT.OaT)  PHODJCEO  hT  RUNNING  ^OhmP 


RECOROET  ON  21-HAY-76  AT  10(27  INTCGN  IN(  « RUN  CPU  TJME 


SLK  NUMBER 
TIME 

0. 0000E<01 
9.0000E-01 
1.0000E*0d 
l.S000E*0v 
2.0000E*0O 
2,S000E*0t> 
3.0000E*0^ 
3.3000E*0i) 
4.3000E*00 
4.5000L*0U 
^.0000E*09 
:>.$300E*0v 
6 • 0000E  * 0 J 
6.5300E*0J 
7,030.'E*04- 
7,5P03E*0* 
8.0'*00E*0«’ 
a.^p03e«0,7 

9,‘>300E*0iJ 
i.0;*00E*0i 
l.0502E*0l 
i.i.’oaf:*?! 
i,i^0eE*0i 
l,2»^00t*0l 
l.2^00E*01 
1.300JE*01 
l,3500E*0l 
l,4P03£*0l 
1.4500E*01 
l,500eE*01 
1 .5500E*01 
i.6K)00E*01 
1.6500E*01 
1 , 7000e*0i 
1.7500E*01 
1.8?00E*01 

1. B500E«01 
1.9000E«01 
1.9500E*01 
2.0000E*01 


« 9 

Y 

f T 

•1.0000E*0I 

-7,802JE*3J 

-3.26V7E*0.’ 

1.0834E*0,T 

4,5316C*P3 
6,32«ftE*03 
6.05b9£*0.T 
4,2;i0E*33 
1.0S/3E*03 
-2.8736E-01 
-2.216&G*03 
-3.2077f *03 
-3.6830E*03 
-3,4467E*3J 
-2,69Xie*03 
-1.7U5£*03 
-7,4723E-ei 
l,3925E-0l 
6,O039E-0l 

1.2940E^00 

1,5817E*00 

1,6747E*03 

1.6303E*03 

1,39^)3E*03 

l,l?21i*03 

7,6292C-ei 

4,1668E-01 

9»5822E-02 

-1.7541C-01 

•3.ei62E-01 

-5,1606E-0l 

-5,7966E-01 

-5,79b6E-0l 

•^♦2735E-01 

•4,3716E-01 

•3,2402E-01 

-2,0227E-01 

•8,4426E-0? 

1,9623E*02 

i.030ie*0i 

1.6200E-01 


4 

Y o-iur  ooT 
^ T/SEr/StC 
2.0203E*^1 
9.1301E*00 
-I . 4i99E*^0 
-3.5>A87E*23 
-6.4ii7E*00 
-8 . 7314E*00 
-6.3973E*00 
-l.8377E*00 
1.1338E*00 
l.67aAE*00 
2.01l8E*r0 

2.9?9S£*00 
2.9307E*00 
2.1354E*00 
8.9477E-01 
-1.0911E-01 
-4.3939E-01 
•6.6944E-01 
-7.9659E-31 
-8.298lE-^l 
-7.8234E-31 
-6.7313E-01 
-^.2314E-ifl 
-3.5334E-21 
-1 .8258E-<;i 
-2.8393E-02 
1.0382E-31 
2.0134E-31 

2.6253E-3X 

2.B924E-31 
2.B807E-2X 
2.5402E-21 
2.0921E-31 
1.5207E-0X 
9.1624E-^2 
3.39X2E-02 
-1.6379E-02 
-5.6362E-02 
-8 . 349«E-i*2 

-V. 844, '•£-32 

-X.010U-0X 


48 

Y UD1 
FT/SLC 

0. 3003t-UX 
7.78B8L«kl0 
9. 3587l*«J0 
7 . 8866L*M3 
3.S216L«k}0 
X.5770t*30 

-2.4443L*0*1 
-4.93&6t*U0 
~4.6169t*03 
-3 . 9097l*00 
-2.9791L*li)0 
-1 . 7008L*00 
-2,0269L-31 
X.3769L*00 
X.e364L«30 
X . 9843L«00 
X.6432L*d0 
X.S6X6L»00 

1. X911L*00 
7. 9047L-01 
3 . 7423L-31 
8.01'  *5t-03 

-2.9?4it.*ji 

-5.XlN9t-01 
-6. 4564L-31 
-6 . 9694L-31 
-6 . 7622t-3l 
•P.  9843t»*'l 
-4 , 8392L-tfl 
-3.41481-01 
-1  . 9663L-31 
-6.03b9L-02 
5 . 65'>2L-02 
1.47l7t-01 
2 . 3810L-01 
2 . 3925L-ri 
2.43221-01 
2.24!>8L-01 
1.09X2L-01 
X.4309L-01 
9.2^381-02 


la 

8L<t  10 

-i.0a0?E*02 
•9.1228£*01 
-1.1619E*01 

2 . X669£*00 

2. X3X6E*01 
4.0603E*01 
3.6825E*01 
1.8ll?E*01 

3.  /X46£*0a 
-P. 7472E-0X 
-4,X30fiE*00 
-X,1248E*01 
-X.439eE*01 
-1.2683£*01 
-8 . X466E*00 
-3 . 423?E*00 
-1 .4944E*00 

2.1851E-01 
X .9308E*00 

2. P98lE*00 
3.1633£*00 
3.3498£«00 
3.230'>£«00 
2 . 79P7E*00 
2.2.142E*00 
l.P258E*00 

H.  J337E-01 

I .  Y164E-01 
-3.5382E-01 
-7.ft323E-01 
•1 .0321E*00 
-1.1593E*00 
-1.159lE*00 
-X.0947£*00 
-0, /433E-01 
-b.4905E-0l 
-4.0484£-01 

-1 .6985E-01 

3. V241E-02 

2.0602E-01 

3.2399E-01 


1,36  StC. 
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APPKNDIX  C — continued 

(2)  (Jraphica!  Output  {File  OVTFT.PPL  printed  using  'P'  switch) 
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APPF.NDIX  1) 

Main  Program  F.xaniple  for  Creating  an  Input  File  for  TRANS  without  Running  BOMMP 


C MAIN  PROGRAM 

EXAMPLE  OF  A MAIN  PROGRAM  WHICH  USEE  A PACKAGE  OF 
SUBPROGRAMS  THAT  WRITE  OUTPUT  INTO  A BINARY  FILE  SUITABLE 
FOR  PROCESSING  BY  THE  PROGRAM  •’TRANS" 

THE  SUBPROGRAM  PACKAGE  CONSISTS  OF  THE  FORTRAN 
FILE  PAKAGE.Fa  AND  THE  MACRO  FILE  CPU. MAC 

THE  fortran  file  CONTAINS  I 

SUBROUTINE  OPI N I T ( COMENT ) 

- RECORDS  parameters  AND  LABELLING  INFORMATION 

SUBROUTINE  OUTPAKtCLAST.C) 

- RECORDS  OUTPUT  FOR  EACH  INDEPENDENT  VARIABLE  STEP 

FUNCTION  PAKBlNt VALLE) 

- USED  BY  OUTPAK 

SUBROUTINE  F I NB I N< TCPU. LOU ) 

- RECORDS  CPU  TIME  AND  FINALIZES  THE  OUTPUT  FILE 

THE  macro  file  CONTAINS  ; 

SUBROUTINE  CPUt TCPU) 

flRSl  CALL  PROVIDES  CPU  TIME  SINCE  JOB  LOGGED  IN» 
SUBSEOUENT  CALLS  PROVIDE  CPU  TIME  SINCE  PREVIOUS 
CALL 

ARRAYS  WITH  VARIABLE  DIMENSIONS  AREt  CLASTtNBLK). 
C(NPLK)#  AND  COMENTt 2»NBLK > , WHERE  NBLK  IS  THE  NUMBER  OF 
OUTPUT  VARIABLES-  AT  EACH  INDEPENDANT  VARIABLE  VALUE.  PRIOR 
TO  RECORDING  THE  OUTPUT  VALUES.  THE  INDEPENDANT  VARIABLE 
MUST  BE  STORED  IN  "T"  AND  THE  OUTPUT  VALUES  MUST  BE  STORED 
IN  THE  ARRAY  "C".  THE  ARRAY  "CLAST"  IS  USED  INTERNALLY  BY 
THE  SUBROUTINE  "OUTPAK"  TO  STORE  THE  LAST  RECORDED  OUTPUT 
VALUES  AND  THUS  SHOULD  NOT  BE  SET  BY  THE  USER 

FOR  THIS  EXAMPLE.  WHERE  NBLK  = 7.  THE  FOLLOWING 
STATEMENTS  MUST  BE  INCLUDED 

DIMENSION  COMENTt 1a).CLAST(7).C<7).NAM(2> 

COMMON  /ZOUTZ/LOU.NBUFF. BUFFERt  I 26 ) . NFLAG. X PC. NBLK 
COMMON  /XOUTX/TI TLEC 12) . TZERO. TSAMP. TTOT.  DT.  T.XLABtP) 

THE  FOLLOWING  VARIABLES  ARE  SET  IN  THE  DATA 
STATEMENTS  BELOW: 

XPC  - PARAMETER  SUCH  THAT  A CURRENT  OUTPUT 

VALUE  (Yl)  IS  RECORDED  ONLY  IF  ABS(YI-Y0) 

IS  GREATER  THAN  ( 0 . 0 1 • XPC* Y 0) . WHERE  Y0 
IS  THE  LAST  VALUE  OF  THE  OUTPUT  VARIABLE 
STORED  AT  A PREVIOUS  TIME  STEP 
NBLK  - NUMBER  OF  OUTPUT  VARIABLES 
LOU  - LOGICAL  UNIT  NUMBER  FOR  THE  OUTPUT  FILE 
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APPENDIX  I) — continued 


NAMte)  - OUTPUT  FILENAME  (MUST  HAVE  EXTENSION 
NAME  "DAT") 

TITUEdP)  - TITLE  OF  UP  TO  60  CHARACTERS 
COMENTCNBLK)  - LABELS  OF  UP  TO  10  CHARACTERS  PER 
OUTPUT  variable.  THE  ARRAY  IS  SET  UP  SUCH 
THAT  FOR  VARIABLE  K.  COMENT( 2*K- I ) AND 
C0MENT(2*K>  CONTAIN  THE  FIRST  AND  SECOND 
GROUPS  OF  FIVE  CHARACTERS  RESPECTIVELY 
DT  - REFERED  TO  IN  "TRANS"  AS  THE  INTEGRATION 

INTERVALl  ITS  ONLY  PURPOSE  IS  TO  IDENTIFY 
THE  DATA 

TZERO  - LOWER  INDEPENDANT  VARIABLE  LIMIT 
TSAMP  - INDEPENDANT  VARIABLE  INTERVAL 

TTOT  - UPPER  INDEPENDANT  VARIABLE  LIMIT 

NFLAG  - A FLAG  SET  EQUAL  TO  1 INITIALLY 
XLAB(2)  " INDEPENDENT  VARIABLE  LABEL  OF  UP  TO  10 
CHARACTERS 

DATA  XPC.NBLK. LOU/0. ei>7. 1/ 

DATA  (NAM(J)» J*l,2)/ ‘FUNCT. DATS'/ 

DATA  (TITLE(J).’j«1.9>/'  VARIOUS  FUNCTIONS 

I OF  X'/ 

DATA  (COMENT(J>,J=l, 14)/ • EXP<2X)  SDUARE 

I*  SIN(8X>  '.'XwSINtSX)  •SIN(8X**2) • SAW 
2 • NO  RMAL  ■ / 

DATA  DT. TZERO. TSAMP, TTOT/ 0.0. 0. 0, 0 . 004, 2 . 0/ 

DATA  NFLAG/ I/, (XLAB( J), J* 1. 2)/ • X '/ 

OPEN  0/P  FILE  I RECORD  THE  LABELLING  INFORMATION 

OPEN(UNI T»lOU. FILE=NAM) 

CALL  OPINIT(COMENT) 

SUBROUTINE  CPU(TCPU)  - FIRST  CALL 
CALL  CPU(TCPU) 

TCPU  contains  CPU  TIME  SINCE  LOGGED  IN 

INITIALIZE  INDEPENDENT  VARIABLE 
T»TZERO 


TIME  LOOP  IN  WHICH  OUTPUT  VALUES  ARE  CALCULATED 
AND  THEN  RECORDED  IN  THE  OUTPUT  BINARY  FILE  WITH  THE 
USE  OF  THE  SUBPROGRAM  PACKAGE 

CONTINUE 

X«T 


FOR  THE  INDEPENDENT  VARIABLE  X.  THE  FOLLOWING 
7 FUNCTIONS  ARE  CALCULATED  AND  RECORDED 

C (1)  EXP<2X) 

C< I )»EXP(2.0*X) 


APPENDIX  D — continued 


C ce>  SflUARE  MAVC  WITH  PERIOD  *.4 

PERIOD-D.4 

ircx.GT.TZCRO>  eo  TO  IID 

cte>>i.D 

SMI TCH-TSAHP>PERIOO/e. I E>« 

110  irtSWITCH.LT.e)  GO  TO  140 
lF<C<e).GT.0.0)  GO  TO  ISO 

c<e)«i .0 

GO  TO  130 
120  C<2)a>|.0 

130  SWITCHsSHITCM-PERlOD/2.0 
GO  TO  MO 

140  SMITCM«SUITCH«TSAMP 

C C3)  SIN(8X> 

C( 3>>SIN(8.04X> 

C <4>  X*SIN(8X) 

C( 4)aX«CC  3) 

C (5>  SIN(8X**2> 

C( S)-S1N(8.0*X**2) 

C (6)  SAW  TOOTH  WITH  PERIOD  0.4 

IFCX.GT.TZERO)  GO  TO  I S0 

c(6>>e.e 

GO  TO  170 
150  C(6>>X 

160  1F<C(6).LT. PERIOD)  C«  TO  170 

C<6>«C(6>-PERI00 
GO  TO  160 
170  CONTINUE 

C <7>  STANDARD  NORMAL  CURVE 

C( 7)a0. 3989* EXP( -X** 2/2.0) 

C OUTPUT  VARIABLES  AT  EACH  ITERATION 

CALL  OUTPAKtCLAST.C) 

T»T*TSAMP 

IFtT.LE.TTOT»0.5*TSAMP)  GO  TO  10 


C OBTAIN  CPU  TINE  SINCE  LAST  CALL  COF  CPU)*  OUTPUT  LAST 

C BUFFERFULL*  AND  CLOSE  OUTPUT  FILE 

CALL  CPU(TCPU) 

CALL  FINBIN(TCPU«LOU) 

CLOSE(UNIT>LOU.FILE>NAM) 

END 
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APPENDIX  t: 

Parameters  and  l.abellin)>  Information  at  Bet’inninj’  of  the  Input  File 


Word  Niiniber(s) 

I 

2 13 
14  17 

IS 

14 

20 

21 

22  (24  3N) 


Contents 


Number  of  output  values  (N) 

Title  of  60  characters 

Pate  and  time  of  creation  of  the  output  file  obtained  when  running 
BOMMP.  which  is  the  input  tile  for  TRANS 

Output  timet  interval 

Integration  interval  used  when  running  BOMMP  (for  reference  only) 

I’pper  time+  value 

1 ower  timet  salue 

N I groups  ot  three  words,  where  the  lirst  word  contains  a bliKk 
number  and  the  second  and  third  words  contain  the  corresponding 
block  label  ol  ten  characters  specified  when  running  BOMMP 
(blank  il  not  specified).  The  first  group  of  three  words  contains  the 
indc|X'ndent  variable  blix'k  number  ( — 1)  and  label 


t or  independent  variable. 


I 

! 


APPENDIX  F 

(/)  Output  Proftnim  Structure:  Mujor  Subroutines 

Names 

Ellcct 

Section(s)  for 

Reference 

TRANl 

Initiates  control  of  the  projiram 

2 

TRAN2 

Controls  the  How  of  the  program 

2 

UNCOM 

Unpacks  block  labels 

5.1 

INITIA 

Unpacks  initial  set  of  block  values 

5.1 

UNPACK 

Unpacks  subsequent  sets  of  block  values 

5.1 

PNTHDS 

Produces  headings  for  Telet>pe  and  line  printer 
tabular  output 

2.2.  3.4 

PRINT 

Produces  Teletype  and  line  printer  tabular  output 

2.2.  3.4 

HEADNS 

Produces  headings  for  Teletype  and  line  printer 
graphical  output 

2.2.  3,4 

I PI  o r 

Produces  Teletype  and  line  printer  graphical  output 

2.2.  3.4 

CAXIS 

Calculates  origins  and  scales  for  incremental  plotter 
output 

3.5 

STRIP 

Produces  curves  for  incremental  plotter  output 

3.5 

GRID.  PLTHD 

Produces  axes  and  labelling  information  for  incre- 
mental plotter  output 

.3.5 

Al'l’KNDIX  I-  -continui'tl 

(-)  Oiilpiil  I’rdi’nim  Slrm  liird  AVn /n-  Suhri>i4lint  \ iiiut  / i-'ii  luins 


Nanie(s) 

1 Iloa 

Section(s)  for 
Reference 

DlNCi 

SouikIn  lu-ll  I'll  the  lclcl\|X- 

t XPANDi 

1 xpaiuls  ,in  .irra\  in  the  ilclmcil  number  nt  dcmcnt> 
(x'c  ReleiciKe  2) 

I cure  K i 

Searches  for  a lilc  on  the  disk 

GTAG,  INBl  K. 

ISRT 

Accepts  lists  of  block  numbers  Irom  the  Eelelcpe  for 
each  form  of  output 

ISRC  H 

Determines  position  of  block  m the  list 

MAXMIN,  XSt  Al  1 , 

^ SCALF;.  (MIN 

Calculates  coiuemenl  losver  aiul  upper  limits  for 
graphical  output 

.7  3 

MXM.  UPDATI 

Accepts  new  sc.tle  limits  from  the  Eelet\pe  and  updates 
the  file  1 R \ \S  SC  1 

3.2 

NBR 

Used  bv  INlI  fA'  and  'UNPXC  K'  to  read  in  a new 
record  of  I2('  dat.i  words 

.3.1 

NUMBE  R.  SVMBOI  + 

I’lols  numbers  and  semhoK  on  the  incremental  plotter 

3.5 

ovrRj.A,  OUST 

Used  tor  'o\crl;i\'  plots 

3.5 

in.oT* 

I’lots  lines  and  points  on  the  incremental  plotter 

3.5 

RIX  1 

Used  to  tielermine  whether  a tile  m.i\  he  read  from  the 
llisk 

RITCPU 

Outputs  CPU  run  time  on  line  printer  output 

3.4 

TI7 

E nables  either  Eeletype  or  line  printer  output  to  be 
specified 

2.2 

TIMOLIT 

Outputs  C'E’U  run  titiie  to  Teletype 

2.1 

TXT.  APE  Nl). 

OB  PAIN 

Accepts  labelling;  modifications  from  the  Teletype  and 
updates  the  file  E RANS.I  Bl. 

3.2 

XAXIS,  YAXIS, 

C'OMC  EEK.  NIC  KI, 
NTCK2 

Used  by  "CiRID'  to  plot  axes  and  labelling  information 

3,5 

ZE  RO,  RE  MVE; 

Removes  block  numbers  from  lists  when  output  is 
completed 

* ARL  I ibnirv  subroulinc. 
^ Writlcn  in  MAC  RO-IO. 
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Kf  f,  . 

^ r 

r T 

0 Y DC’’ 

F T/Sf L 

^ f OPlF  DO’ 
f VCf.L^CF.t 


1 IMF. 
jfX 


NCN-lINERP  CPFilNC,  HUM  D»liHPGT 
DfiTp  f^TCRFD  IN  OUTP’.nfJT)  PRODUCF.D  Pr  PUNNING  POWP 


•Mo.j  7p  i0.^^9  ^ \ .mm\:  01  cpu  timu  - 1.3c  cfx. 


FIG.  1b.  INCREMENTAL  PLOTTER  OUTPUT;'OVERLAY'plOTS 
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